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(Dated: February 7, 2008) 

Abstract 

Using 29.1 fb^ 1 of data accumulated at the T(45) with the Belle detector at KEKB, we have 
studied the decay modes B° — > A+p7r + 7r~, B~ — > Afpir~~, and B° — > Afp. We report branching 
fractions of exclusive B decays to charmed baryons with four-, three- and two-body final states, 
including intermediate S++ and £° states. We observed B° — > S c (2455) ++ p7r~ for the first time 
with a branching fraction of (2.38I i Iq;|| ±0.41 ±0.62) x 10~ 4 and observed evidence for the two-body 
decay B~ -> £ c (2455)°p with a branching fraction of (0.45l[j;?g ± 0.07 ± 0.12) x 10" 4 . We also set 
improved upper limits for the two-body decays B° — ► Afp and B~ — > E c (2520)°p. 

PACS numbers: 13.25.Hw, 14.20.Lq 
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Baryon production in flavored meson decays is unique to the B meson system due to 
the heavy mass of the constituent b-quark. Several studies of inclusive charmed baryon 
production in B meson decays JIJ have been made and a large branching fraction for B — > 
A^X of (6.4±1.1) % has been reported. However, the mechanism is not well understood. The 
measured inclusive A+ momentum spectra indicate that multi-body final states are dominant 
in baryonic B decays. With a data sample of 2.39 fb -1 , CLEO has studied exclusive 
charmed baryonic decay modes and measured the branching fractions for B° — > A^p7i + n^ 
and B~ — > A+p7r~. They found no evidence for B° — ► A^p and provided an upper limit. So 
far, no observations of two-body decays have been reported. On the other hand, there are 
theoretical predictions for branching fractions of two-body baryonic modes based on a pole 
model |§, a QCD sum rule [|J, a diquark model ||, and a bag model ||. The predictions of 
the different models vary by an order of magnitude, and experimental measurement can be 
used to discriminate among them. We have made a systematic study of exclusive charmed 
baryonic decays of B° and B~ mesons into four-, three- and two-body final states including 
intermediate resonances, by analyzing the A^pn + n^, A+p7r~ and A+p final states. 
Charge conjugate modes are included unless otherwise mentioned. This analysis is based on 
a data sample of 29.1 fb -1 corresponding to 3.17x 10 7 BB pairs. The data were accumulated 
at the T(4S) resonance with the Belle detector at the KEKB asymmetric collider of 3.5 GeV 
e+ and 8.0 GeV e~ 0. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a three- 
layer silicon vertex detector (SVD), a 50-layer cylindrical drift chamber (CDC), a mosaic of 
aerogel threshold Cerenkov counters (ACC), a barrel-like array of time-of-flight scintillation 
counters (TOF), and an array of CsI(Tl) crystals (ECL) located inside a superconducting 
solenoid coil that provides a 1.5 T magnetic field. An iron flux return located outside the 
coil is instrumented to detect muons and Kl mesons (KLM). The detector is described in 
detail elsewhere ||. We use a GEANT based Monte Carlo (MC) simulation to model the 
response of the detector and determine the acceptance ||. 

In searches for the decay modes B° — > A+p 7r + 7r~, B~ — > A+p tt~ , and B° — > A+p, 
the A+ — > pK~7i + decay mode is used. Particle identification information from the CDC 
dE/dx, ACC and TOF is used to provide a mass assignment for each track. A likelihood 
ratio LR(A,B) = La/ (La + Lb) > 0.6 is required to identify a particle as type A, where 
B is the other possible assignment among n^, and p(p). Electron and muon candidate 
tracks are removed if their probabilities from the ECL, CDC dE/dx and KLM are greater 
than 95%. Candidate A+'s are tagged if the invariant mass of the p, K~ and tt + track 
combination is within 0.010 GeV/c 2 of the A+ mass; tagged events are then examined for 
the three search modes by adding p, ir~, and tt + tracks. The width a A + is found to be 
4.9 MeV/c 2 , consistent with the MC. 

In order to select B meson candidates, we use the beam energy-constrained mass and 

energy difference, which are defined as Mb c = \J El c&m — (J2pl) 2 an d AE = Ei ~ E^ eam 
in the center-of-mass (CM) frame of the e + e~ collision. -Ebeam is the beam energy, and £7j 
and pi are the energy and momentum vector for the i-th daughter particle of a B candidate. 
B candidates are selected with a loose cut to retain sideband events by requiring Mb c > 
5.2GeV/c 2 and \AE\ < 0.2 GeV. A vertex-constrained fit for the three daughter tracks is 
carried out at the A+ vertex. For each decay mode, the virtual A+ track and additional 
tracks are required to form a good vertex. If there are multiple candidates for both Ajt and 
B, the candidate with the minimum \ 2 — X 2 A + + Xb + (^bc — 5.279) 2 /cx| fbc is selected. Here, 
X K + and Xb are ^ ne X 2 ' s from the fits for the A+ and B vertices, respectively, and cr Mhc is 
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the MC value of the M DC width (2.8 MeV/c 2 ). Loose cuts on \\+ an d Xb are applied to 
remove background from tracks arising from Kg and A decays. 

Event selection requirements are optimized using signal MC events and continuum back- 
ground MC events consisting of uu, dd, ss, and cc quark-antiquark pairs generated with the 
expected fractions. To suppress the continuum background, we use a Fisher discriminant 
constructed from 10 variables: 8 modified Fox Wolfram moments | 10| , cosO#, and cos0 A +. 



Here, cosO# is the cosine of the direction of the B meson with respect to the electron beam 
direction, and cos9 A + is the cosine of the direction of the daughter A+ with respect to the 
thrust axis of the tracks not associated with the B candidates. Both quantities are defined 
in the CM system. A set of 10 coefficients for each mode is optimized to maximize separa- 
tion of the signal from the continuum background. The probability density functions for the 
signals and for the continuum, P slg and P con respectively, are parameterized with Gaussian 
functions for the three search modes and for the continuum events. A cut on the likelihood 
ratio _R s f w = -P S ig/(-P S ig + -Peon) > 0.6 is applied to all decay modes. In the MC simulation 
this cut removed 76% of the continuum background while retaining 86% of the signal for 

A+p7T + 7T~. 

Figure |I| shows the M DC and AE distributions for the three decay modes, after a tight 
cut is made in the (AE,M\, C ) variable not plotted. The M DC background distributions are 
parameterized by the ARGUS function while a Gaussian is used for the signal. The 
AE distributions are fitted with a second-order polynomial for the background and a double 
Gaussian for the signal. Here, the width parameters are fixed to the values fitted to the 
signal MC events. The mean and width of M DC in the data are found to be consistent with 
the MC values of 5.279 GeV/c 2 and 2.8 MeV/c 2 , respectively. The width of AE is also 
consistent with the MC value (9.9 MeV) when fit to a single Gaussian. We obtain signal 
yields of 154ljg and 38. 8^7^ from the fits to the M^ c distributions (a) and (c), and 1411^ 
and 30.2ig 4 from the fits to the AE distributions (b) and (d), respectively. Here, we choose 
the asymmetric range of —0.100 < AE < 0.200 GeV to exclude feed-down from higher 
multiplicity modes with extra pions; these produce the structure observed in the region 
AE < —0.150 GeV. Since M DC is used in the \ 2 calculation for the best candidate selection 
as described previously, we use the yields resulting from the fits to the AE distributions to 
calculate branching fractions. 

We observe B° — > A+j>7r + 7r~ and B~ — > A^pn~ signals. For B° — > A+p we find a 
statistical significance of only 1.9 a from a fit to a Gaussian function for the signal with 
mean and width fixed to those from the signal MC, and a linear background function. We 
thus set an upper limit of 6.1 events at the 90% confidence level based on the likelihood 
function, using the Bayesian method with a prior uniform in the branching fraction. 

Table | summarizes the observed yields and branching fractions. Here, the detection 
efficiencies are calculated assuming nonresonant decays and do not include the branching 
fraction £>(A+ — > pK~7r + ) = (5.0 ± 1.3)% [□[]. We assume the fractions of charged and neu- 
tral B mesons to be equal in the branching fraction calculations. We include a correlated 
systematic error of 2% per track for tracking and particle identification. Systematics due to 
the Xa + j Xb an d -Rsfw cuts are estimated by varying cut values. The signal shape systematic 
error is evaluated from the variation in fit results obtained with different-order polynomials 
used for the background and single and double Gaussians used for the signal. The resulting 
total systematic errors for A+p7i + 7i~ , A+pTi~ and lytp are 17.2%, 14.8% and 13.3%, respec- 
tively. Table | shows the CLEO measurements renormalized to the same B(A+ — > pK~7T + ) 
for comparison. Our branching fraction for B° — > A+p7r + 7r~ is consistent with their mea- 
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FIG. 1: M hc distributions for \AE\ < 0.030 GeV and AE distributions for M bc > 5.270 GeV/c 2 : 
(a) and (b) for B° —> A+p7r + 7r~, (c) and (d) for B~ — > A+pir~, and (e) and (f) for B° — > A+p. 
Points with errors indicate the data and the curves indicate fits (see text for details). 



surement; however, our result for B — ► A + p7r 
more restrictive upper limit on B° — > A + p. 



is somewhat lower (1.5 a). We also set a 



TABLE I: Branching fractions for B° — > A + p7r + 7r~, B~ — > A + p-7r _ , and -B — ► A + p. The errors are 
statistical, systematic, and a common error due to the uncertainty in the value of B(A+ — > pK~ir + ). 
The CLEO results are renormalized to $(A + — > pK~ir + ) = (5.0 ± 1.3)% [|l2| for comparison. 



Mode 



Efficiency 



Yield Significance 



B fxlO" 4 ) 



CLEO fxlO" 4 ) 



B° - 


-> A + p7T + 7T 


8.07 


141+H 


12.2 


11.0±1;| ± 1.9 ±2.9 


B- - 


-> A+p7T~ 


10.2 


30-2+™ 


6.0 


1.8T±g;|g ± 0.28 ± 0.49 


B° - 




12.9 


2 4+ 2 ' 1 


1.9 


0.12±g;J? ± 0.02 ± 0.03 



7+4.0 
'-3.7 



±2.7 ±3.0 



< 6.1 (90% CL) 



< 0.31 (90% CL) < 1.85 (90% CL) 



Figure fj shows the A^n^ 1 invariant mass distributions in the B signal region, \AE\ < 
0.030 GeV and M bc > 5.270 GeV/c 2 . Significant signals are observed for the £ c (2455) 
and E c (2520). The shaded histograms are the distributions for events in the sideband region 
0.040 < \AE\ < 0.100 GeV, normalized to the signal region \ AE\ < 0.030 GeV; these account 
for continuum S c background. The two curves indicate the results of separate fits to the 
distributions for the B signal and the sideband regions, with S c masses and widths fixed to 
fit values for the signal MC events generated with PDG values for masses and widths |12| . 
The background shapes are taken from a nonresonant signal MC. To extract the E c yields, 
we performed a simultaneous likelihood fit to the distributions for the B signal and sideband 
regions. We express the expected number N% c of B events as iV 2c = N Bb — r ■ N sb , where 
Nsb is the yield in the B signal region, N s b is the yield in the sideband region, and r = 0.5 is 
the normalization factor due to the ratio of their AE ranges, assuming a linear background 
shape. 
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FIG. 2: Invariant mass distributions (a) M(A+7r + ) and (b) M(A+7r~) for B° — > A+p7r + 7r~, 
and (c) M(A+7r~) for B — > A+p7r~. Points with errors and shaded histograms indicate the 
distributions for the B signal and the sideband regions, respectively. The curves indicate fits (see 
text for details). 

Table summarizes the observed signal yields and branching fractions. We observe the 



TABLE II: Efficiencies, yields, significances and branching fractions for decay modes with 
resonances. The errors are statistical, systematic, and a common error due to the uncertainty in 
the value of #(A+ -> pK~ir + ). 



Mode 




Efficiency (%) 




Yield Si; 


mificance 


B (xlO- 4 ) 


B ^ 


£ c (2455) ++ p7T- 


4.93 




18.6+tl 


5.3 


2.38±g;f| ± 0.41 ± 0.62 


B ^ 


£ c (2520) ++ p7T- 


6.38 




16.5^:1 


3.5 


1.63±g;|J ± 0.28 ± 0.42 


B ^ 


S c (2455)°p7r+ 


4.80 


< 


6.4i|;? 
11.6 (90% CL) 


2.6 


0.84±g;|| ± 0.14 ±0.22 
< 1.59 (90% CL) 


B ^ 


£ c (2520)Vr+ 


6.35 


< 


A S+ 4 - 5 
4 -8_ 4 . 

11.7 (90% CL) 


1.2 


0.48±g;^ ± 0.08 ± 0.12 
< 1.21 (90% CL) 


B~ - 


£ c (2455)°p 


6.00 


< 


^-1.8 

8.5 (90% CL) 


3.0 


0.45±g;53 ± 0.07 ± 0.12 
< 0.93 (90% CL) 


B- - 


£ c (2520)°p 


7.47 


< 


i 7 +i.8 
5.2 (90% CL) 


1.8 


0.14^;^ ±0.02 ±0.04 
< 0.46 (90% CL) 



B° — ► £ c (2455) ++ p7r~ decay for the first time with a statistical significance of 5.3 a. We also 
see 3.5a evidence for B° — > £ c (2520) ++ ]97r~, 2.6 a evidence for B° — > E c (2455)°p7r + , and 
less evidence for B° — > E c (2520)°p 7r+ - We see 3.0<r evidence for the two-body decay B~ — > 
S c (2455)°p, and less evidence for B — > S c (2520)°p. For those modes with a significance of 
three sigmas or less, we set upper limits on their branching fractions. 
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Our results provide stringent constraints upon theoretical predictions 0, f|, [5], |[. The 
predictions for B° — > A+p in |3], [|, [5| were already much larger than the CLEO experimental 
upper limit |2[]; here we set an even more restrictive upper limit. A recent study based 
on a bag model || gives predictions of branching fractions of < (0.1 ~ 0.3) x 10~ 4 for 
B° — > A+p and (4.3 ~ 15.1) x 10 -4 for B~ — > A+p7r~. Our upper limit for B° — > A+p does 
not contradict this model, while our measured result for B~ — > A+p7r~ is much smaller than 
its predicted value. 

In summary, we have observed the exclusive three-body decay B° — > £ c (2455) ++ ]97r~ for 
the first time and observed evidence for the exclusive two-body decay B~ — > S c (2455)°p. 
We make improved measurements of the branching fractions for B° — > A+]97r + 7r~ and B~ — ► 
A+p7r~, and also set a more restrictive upper limit on B° — ► A+p. 
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